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' Abstract 



In this paper the asymptotic behavior of an unstable integer-valued autoregressive 
model of order p (INAR(p)) is described. Under a natural assumption it is proved that 
the sequence of appropriately scaled random step functions formed from an unstable 
INAR(p) process converges weakly towards a squared Bessel process. We note that this 
■ limit behavior is quite different from that of familiar unstable autoregressive processes of 

c| . order p. 

1 Introduction 

>. 

Recently, there has been remarkable interest in integer-valued time series models and a num- 
ber of results are now available in specialized monographs (e.g., MacDonald and Zucchini [32], 
Cameron and Trivedi [IT] , and Steutel and van Harn [58]) and review papers (e.g., McKenzie 
6J, Jung and Tremayne [33], and WeiB [60J). Reasons to introduce discrete data models come 



OO 

o . 

ON . from the need to account for the discrete nature of certain data sets, often counts of events, 
objects or individuals. Examples of applications can be found in the analysis of time series of 



count data on the area of financial mathematics by analyzing stock transactions (Quoreshi |51j). 
insurance by modeling claim counts (Gourieroux and Jasiak [23]), medicine by investigating 
disease incidence (Cardinal et al. [12]), neurobiology by change-point analysis of neuron spike 
train data (Belisle et al. [3]), optimal alarm systems (Monteiro et al. [31]), psychometrics by 
treating longitudinal count data (Bockenholt [6], [7]), environmetrics by analyzing rainfall mea- 
surements (Thyregod et al. [59]), experimental biology (Zhou and Basawa [61]), and queueing 
systems (Ahn et al. [1] and Pickands III and Stine [50]). 

Among the most successful integer-valued time series models proposed in the literature 
we mention the INteger- valued AutoRegressive model of order p (INAR(p)). This model 
was first introduced by McKenzie [35] and Al-Osh and Alzaid [2] for the case p — 1. The 
INAR(l) model has been investigated by several authors. Franke and Seligmann [19] analyzed 
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maximum likelihood estimation of parameters under Poisson innovation. Du and Li jTTJ and 
Freeland and McCabe [2T] derived the limit-distribution of the ordinary least squares estimator 
of the autoregressive parameter. Silva and Oliveira [51] proposed a frequency domain based 
estimator, Brannas and Hellstrom [8] investigated generalized method of moment estimation, 
Silva and Silva [56] considered a Yule- Walker estimator. Jung et al. [32] analyzed the finite 
sample behavior of several estimators by a Monte Carlo study. Ispany et al. [27], [28] derived 
asymptotic inference for nearly unstable INAR(l) models which has been refined by Drost et 
al. [15] later. A Poisson limit theorem has been proved for an inhomogeneous nearly critical 
INAR(l) model by Gyorfi et al. [21]. 

The more general INAR(p) processes were first introduced by Al-Osh and Alzaid [3] . In their 
setup the autocorrelation structure of the process corresponds to that of an ARMA(p,p — 1) 
process, see also Section [2j Another definition of an INAR(p) process was proposed indepen- 
dently by Du and Li [T7] and by Gauthier and Latour [22| and Latour [3JJ], and is different 
from that of Alzaid and Al-Osh [3J. In Du and Li's setup the autocorrelation structure of an 
INAR(p) process is the same as that of an AR(p) process. The setup of Du and Li [17] has been 
followed by most of the authors, and our approach will also be the same, see Section [2j The 
INAR(p) model has been investigated by several authors from different points of views. Drost 
et al. [2] provided asymptotically efficient estimator for the parameters. Silva and Oliveira [55] 
described the higher order moments and cumulants of INAR(p) processes, and Silva and Silva 
[56] derived asymptotic distribution of the Yule- Walker estimator. Drost et al. [16j considered 
semiparametric INAR(p) models and proposed efficient estimation for the autoregression pa- 
rameters and innovation distributions. Recently, the so called p-order Rounded INteger- valued 
AutoRegressive (RINAR(p)) time series model was introduced and studied by Kachour and Yao 
[33] and Kachour [31]. The broad scope of the empirical literature in which INAR models are 
applied indicates its relevance. Examples of such applications include Franke and Seligmann 
[T9] (epileptic seizure counts), Bockenholt [7J (longitudional count data), Thyregod et al. [59] 
(rainfall measurements), Brannas and Hellstrom [8] and Rudholm |53j (economics), Brannas 
and Shahiduzzaman [U] (finance), Gourieroux and Jasiak [23J (insurance), Pavlopoulos and 
Karlis [H] (environmental studies) and McCabe et al. [H] (finance and mortality). 

An interesting problem, which has not yet been addressed for INAR(p) models, is to in- 
vestigate the asymptotic behavior of unstable INAR(p) processes, i.e., when the characteristic 
polynomial has a unit root. In this paper we give a complete description of this limit behavior. 
In particular, it will turn out that an INAR(p) model is unstable if and only if the sum of 
its autoregression parameters equals 1, and in this case the only unit root of the character- 
istic polynomial is 1 with multiplicity one. For the sake of convenience, we suppose that the 
process starts from zero. Without loss of generality, we may suppose that the pth autoregres- 
sion parameter is strictly positive and that the greatest common divisor of the strictly positive 
autoregression parameters is 1, see Remark 12.21 Under the assumption that the second mo- 
ment of the innovation distribution is finite, we prove that the sequence of appropriately scaled 
random step functions formed from an unstable INAR(p) process converges weakly towards a 
squared Bessel process. This limit process is a continuous branching process also known as 
square-root process or Cox-Ingersoll-Ross process. We should also note that the behavior of 
unstable INAR(p) models is completely different from that of familiar (real-valued) unstable 
AR(p) models in at least two senses. On one hand, the characteristic polynomial of an unstable 
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INAR(p) model has only one unit root, namely 1, with multiplicity one, whereas for a familiar 
AR(p) model it may have real or complex unit roots with various different multiplicities. On 
the other hand, in the case of an unstable INAR(p) model there is a limit process which is a 
squared Bessel process, while in the case of an unstable AR(p) model in general there is no 
limit process, only for appropriately transformed and scaled random step functions, see Chan 
and Wei [13], Jeganathan (31] and van der Meer et al. [431 Theorem 3]. 

We remark that our result can be considered as the first step towards the comprehensive 
theory of nonstationary integer-valued time series and investigation of the unit root problem 
of econometrics in the integer-valued setup. Nonstationary time series have been playing an 
important role in both econometric theory and applications over the last 20 years, and a sub- 
stantial literature has been developed in this field. A detailed set of references is given in 
Phillips and Xiao [49] . We note that Ling and Li [40], [H] considered an unstable ARMA 
model with GARCH errors and an unstable fractionally integrated ARMA model. Concerning 
relevance and practical applications of unstable INAR models we note that empirical studies 
show importance of these kind of models. Brannas and Hellstrom [8] reported an INAR(0.98) 
model for the number of private schools, Rudholm [53] considered INAR(0.98) and INAR(0.99) 
models for the number of Swedish generic-pharmaceutical market. Hellstrom [25J focused on 
the testing of unit root in INAR(l) models and provided small sample distributions for the 
Dickey-Fuller test statistic under the null hypothesis of unit root in an INAR(l) model with 
Poisson distributed innovations. To our knowledge a unit root test for general INAR(p) models 
is not known, and from this point of view studying unstable INARfjo) models is an important 
preliminary task. 

The rest of the paper is organized as follows. Section [2] provides a background description of 
basic theoretical results related with INAR(p) models. In Section [3] we describe the asymptotic 
behavior of unstable INAR(p) processes. Under the assumption that the second moment of 
the innovation distribution is finite, we prove that the sequence of appropriately scaled random 
step functions formed from an unstable INAR(p) process converges weakly towards a squared 
Bessel process, see Theorem 13.11 Section [4] contains a proof of our main Theorem 13.11 For the 
proof, we collect some properties of the first and second moments of (not necessarily unstable) 
INAR(p) processes, we recall a useful functional martingale limit theorem and an appropriate 
version of the continuous mapping theorem, see Lemma 15. 1\ Corollary 15. 1[ Theorem 15.11 and 
Lemma 15.21 in Appendix, respectively. 

2 The INAR(p) model 

Let Z + , N, R, R + and C denote the set of non-negative integers, positive integers, real 
numbers, non-negative real numbers and complex numbers, respectively. For all n £ N, let 
us denote by J n the n x n identity matrix. Every random variable will be defined on a fixed 
probability space (Q,A,P). 

One way to obtain models for integer-valued data is to replace multiplication in the con- 
ventional ARMA models in such a way to ensure the integer discreteness of the process and to 
adopt the terms of self-decomposability and stability for integer-valued time series. 
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2.1 Definition. Let (Ek)km be an independent and identically distributed (i.i.d.) sequence 
of non-negative integer-valued random variables, and let ai,...,a p G [0,1]. An INAR(p) 
time series model with coefficients «i, . . . , a p and innovations (sk)keN is a stochastic process 
(X„) n ^_ p+ i given by 

(2.1) x k = + • • • + Yl + £fc ' ke N > 

j=l 3=1 

where for all k G N and i G {1, . . . ,p}, (£k,i,j)jeN a sequence of i.i.d. Bernoulli random 
variables with mean ol{ such that these sequences are mutually independent and independent 
of the sequence (ek)k&N, and Xq, . . . , X_ p+ i are non-negative integer-valued random 

variables independent of the sequences (Ck,i,j)jeH, k G N, i G {1, . . . ,p}, and (e k ) k m- 

The INAR(p) model (12 .ip can be written in another way using the binomial thinning oper- 
ator ao (due to Steutel and van Harn [57]) which we recall now. Let X be a non- negative 
integer-valued random variable. Let (Cj)jen be a sequence of i.i.d. Bernoulli random variables 
with mean a G [0,1]. We assume that the sequence (£j)jeH is independent of X. The 
non-negative integer- valued random variable a o X is defined by 



a o X : = 




if X>0, 
if X = 0. 



The sequence (£i)jsN is called a counting sequence. The INAR(p) model (12 .ip takes the form 

X k = ai o X fe _i H h a p o X fe _p + e k , k G N. 

Note that the above form of the INAR(p) model is quite analogous with a usual AR(p) process 
(another slight link between them is the similarity of some conditional expectations, see (12. 3p ). 
As we noted in the introduction, this definition of the INAR(p) process was proposed indepen- 
dently by Du and Li [TTJ and by Gauthier and Latour [22] and Latour [39], and is different from 
that of Alzaid and Al-Osh [3], which assumes that the conditional distribution of the vector 
(ctioXf, Oi2oX t , . . . , aipoXt) given X t = x t is multinomial with parameters (a%, 012, ■ ■ ■ , ot p , x t ) 
and is independent of the past history of the process. The two different formulations imply 
different second-order structure for the processes: under the first approach, the INAR(p) has 
the same second-order structure as an AR(p) process, whereas under the second one, it has the 
same one as an XRMA(p,p — 1) process. 

An alternative representation of the INAR(p) process as a p-dimensional INAR(l) process 
was obtained by Franke and Subba Rao [20] and see also Latour [38l formula (2.3)]. Accordingly, 
the INAR(p) process defined in (12. ip can be written as 

X k = AoX k _i + e k , keN, 
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where the p-dimensional random vectors X k , s k and the (p x p)-matrix A are defined by 



(2.2) 



Xk 
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a 3 ■ ■ 
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and for a p-dimensional random vector Y = (¥]_, . . . , Y p ) and for a pxp matrix B = J=1 
with entries satisfying ^ fry ^ 1, z,j = 1, . . . ,p, the matricial binomial thinning operation 
B oY is defined as a p-dimensional random vector whose i-th component, % = 1, . . . ,p, is 
given by 

v 

3=1 

where the counting sequences of all 6^ o Yj, i,j = 1, . . . ,p, are assumed independent of each 
other. 

In what follows for the sake of simplicity we consider a zero start INAR(p) process, that is 
we suppose X = X_i = . . . = X_ p+1 = 0. The general case of nonzero initial values may be 
handled in a similar way, but we renounce to consider it. For nonzero initial values the first 
and second order moments of the sequence (X k ) ke z+ have a more complicated form than in 
Lemma 15.11 Further, for proving a corresponding version of our main result (see Theorem 13.11) 
one needs to apply a more general version of Theorem 15.11 which is also valid for random step 
functions not necessarily starting from 0. 

In the sequel, we always assume that E(ef) < oo. Let us denote the mean and variance of 
£ i by /i £ and erf, respectively. 

For all k e Z+, let us denote by J- k the a- algebra generated by the random variables 
X ,X U ..., X k . (Note that ? Q = {0, Q}, since X = 0.) By lEO) . 



(2.3) 

Consequently, 



E(X fc | J- k -i 

E(x fc ; 



Oi\X k -i H h a p X k -p + \i 



keN. 



This can also be written in the form E(X k 

T 



ai E(X fe _i) + • • • + a p E(X k _ p ) + /i £ , k e N. 

AE(X fe _ 1 ) + /i £ e 1; /c e N, where d 



1,0,0, 



,0,0 



opxl 



Consequently, we have 



which implies 
(2.4) 



E(Xt 



fc-i 
i=o 



E(X fc ) = E(e7X A 



fe-i 

/ie^e^A-'ei, 
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Hence the matrix A plays a crucial role in the description of asymptotic behavior of the 
sequence (X fc ) fc ^_ p+1 . Let g(A) denote the spectral radius of A, i.e., the maximum of the 
modulus of the eigenvalues of A. 

In what follows we collect some known facts about the matrix A. First we recall the 
notions of irreducibility and primitivity of a matrix. A matrix M G W xp is called reducible 
if p — 1 and M = 0, or if p ^ 2 and there exist a permutation matrix P G M. pxp and an 
integer r with 1 ^ r ^ p — 1 such that 



P T MP 



B C 
D 



where B G M rxr , D G R^Mp-^ c e R rx(p-r)^ and q e K (p-r)xr ig ft nuU matrix 
A matrix M G M. pxp is called irreducible if it is not reducible, see, e.g., Horn and Johnson 
[26, Definitions 6.2.21 and 6.2.22]. A matrix M G R^ xp is called primitive if it is irreducible 
and has only one eigenvalue of maximum modulus, see, e.g., Horn and Johnson [261 Definition 
8.5.0]. By Horn and Johnson [2"6"l Theorem 8.5.2], a matrix M G IR!J_ xp is primitive if and only 
if there exists a positive integer k such that all the entries of the matrix M k are positive. 

Let us denote by ip the characteristic polynomial of the matrix A, i.e., 

tp(\) := det(A/ p - A) = X p - ai \ p - 1 a p ^X - a p , A G C. 

2.1 Proposition. For a\, . . . , a p G [0, 1], a p > 0, let us consider the matrix A defined in 
(12.21) . Then the following assertions hold: 

(i) The characteristic polynomial tp has just one positive root, g(A) > 0, the nonnegative 
matrix A is irreducible, g(A) is an eigenvalue of A and 



(2.5) ^a fc ^(A)- fe = l, 



k=l 
P 



k=l 



(2.6) ka ^(A)- k = g(A)- p+1 V '(g(A)). 

Further, 

(2.7) 





// d — 1, then A is primitive, g{A) is an eigenvalue of A, the algebraic and 
geometric multiplicity of g(A) equal 1 and the absolute value of the other eigenvalues 
of A are less than g(A). Corresponding to the eigenvalue g{A) there exists a unique 
vector ua G MP with positive coordinates such that Au^ = g{A)uA and the sum of 
the coordinates of ua is 1, namely, ua takes the form 

r 1 T g(A)- i+1 . 1 

U A= u A)1 ,...,u AyP with u Atl :==p A ,_ k+1 , i = ±,..., p. 

1 J l^k=i Qy^) 
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Further, 

(2.8) g(A)~ n A n -»■ II A := u A v T A , as n ->■ oo ; 

where v a 6R p a unique vector with positive coordinates such that A t va = q(A)va 

r iT 

and w^v^ = 1, namely Va takes the form v a = va x, ■ ■ ■ , f An ^i/i 



' - EL^(Ap £ ate(A) - sM-^)) £ tt ' e(A) ' 

/or 2 = l,...,p. Moreover, there exist positive numbers Ca and rA with < 1 
suc/i that for all n G N 

(2.9) ||o(A)~ n A" - n^H ^ cat a, 

where \\B\\ denotes the operator norm of a matrix B £ M pxp defined by := 
SUpiun-! ||-Bx|| . 

Proof, (i): First we check that has just one positive root, which readily yields that 

g(A) > 0. The function Ah> 1- A~*V(A) = ctiA^ 1 H h a p _iA~ p+1 + a p X~ p is strictly 

decreasing and continuous on (0, oo) with lim(l — A _J V(A)) = oo and lim(l — \~ p (p(X)) = 0, 

A^O Atoo 

thus it takes the value 1 at exactly one positive point, which is the only positive root of (p. 



Now we turn to check that A is irreducible. By Brualdi and Cvetkovic [lOj Definition 8.1.1 
and Theorem 1.2.3], a nonnegative matrix B = (&ij)ij=i r .. p is irreducible provided that its 
digraph (directed graph) D(B) (having p vertices labeled by the numbers 1,2, ... ,p and 
an edge from vertex i to vertex j provided £>y > 0) is strongly connected (that is, for each 
pair i and j of distinct vertices, there is a path from i to j and a path from j to i). 
Now a,p > implies that D(A) contains a cycle 1— > p (p — 1)— >■•••— ^ 2 — )■ 1, hence 
-D(A) is strongly connected. 

Using that A is nonnegative and irreducible, by Horn and Johnson Theorem 8.4.4], we 
have q(A) is an eigenvalue of A and hence 

g(A) p - a 1 g(A) p ~ 1 a p _ 1 g(A) -a p — 0, 

which yields (12. 5p . Since 

<p'(X) = pX^ 1 - (p - l)axA p - 2 a p _i, A e C, 

we have 

</(<?(A)) = pg{AY l oc k Q{Ay~ k - J2(P - k)a k g(A) p - k - 1 

k=l k=l 



J2^ k g(A) p - k - 1 = QiAy-^kauQiA)- 



k=l k=l 

which yields (12.61) . 
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Further, (12. 51) yields that 

if <{ = 1, then l = ^a^(A)- fe <^= ^a k . 

k=l k=l 

This readily implies (12. 7ft . 

(ii): By Brualdi and Cvetkovic [TOl Definition 8.2.1 and Theorem 8.2.7], an irreducible nonneg- 
ative matrix B = (bij)ij=i t ... p is primitive provided that the index of imprimitivity of B (the 
greatest common divisor of the lengths of the cycles of its digraph D(B)) equals 1. Now the 
cycles of D(A) are 1— >•«—>• (i — 1) >■ 2 — >■ 1 for all i — 1, . . . , p such that ctj > 

(not considering rotations). Since such a cycle has length z, we get the index of imprimitivity 
of A is d = 1, which yields that A is primitive. 

The other assertions of (ii) except the uniqueness of ua and va follows by the Frobenius- 
Perron theorem, see, e.g., Horn and Johnson [26, Theorems 8.2.11 and 8.5.1]. The uniqueness 
of ua follows by Horn and Johnson [261 Corollary 8.2.6] using that g(A m ) = g(A) m for all 
m G N. The uniqueness of va can be checked as follows. Using that the irreducibility and 
primitivity of A yields the irreducibility and primitivity of A T (see, e.g., page 507 in Horn 
and Johnson [26]), by Horn and Johnson [261 Theorems 8.2.11, 8.5.1 and Corollary 8.2.6] we 
get g(A T ) = g(A) is an eigenvalue of A T , the algebraic and geometric multiplicity of g(A) 
equal 1, corresponding to the eigenvalue g(A) there exists a unique vector va £ K p with 
positive coordinates such that A t va — g(A)vA and the sum of the coordinates of is 
1. Further, by Horn and Johnson [26, page 501, Problem 1], we also have u\va > 0. Using 
that the geometric multiplicity of g(A T ) = g(A) equals 1, we get va '■= u t\ a ^a is a unique 
vector with positive coordinates such that A t va = g(A)vA and u\va = 1. 

The forms of Ua and va can be checked as follows. Using that they are unique it 
remains to verify that the imposed conditions are satisfied by the given forms. We easily have 
Ua has positive coordinates of which the sum is 1. Further, with the notation Aua = 
[(Au A )i, ■ ■ ■ , [Au A )p] , we get 

s EUuig(A)- l+1 g(A) >A . A ., g(A) 

= g(A)u A>1 , 

where the last but one equality follows by (12. 5p . Similarly, for i — 2, . . . ,p, we get 

g(A)~ i+2 

(Au A )i = u A i-i = -, A s_ k+1 = g(A)u A ,i- 

l^k=i Q\ A ) 

Moreover, we easily have va has positive coordinates and 

p / p 



U >* = rfA)^U(A)) S ( e<Ar+1 E a " (Arl 



^ P P 



g(A)-P+^'(g(A)) ^ 



i=i 
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where the last equality follows by (12. 6ft . With the notation A T v A = [(A T v A )i, • • • , (A 7 "^)?] 
we get for all i — 1, . . . ,p — 1, 



E*-i0(^)~ fc+1 / " 

(A = ai«A,i + v A , i+ i = . k ~ l .: ., — I a 



e(A)-p+V(e(^)) 



i-t 



t=i+l 



ELi 9{A 



,-k+l 



e=i+i 



i-i 



g(A)-P+^'(g(A)) 
= g(A)v Aii . 

Finally, using that Y^=i a kg(A)~ k = 1, we get 
(A T v A ) P = a p v A ,i = a P 



g(A)-P+^'(g(A)) 



i=i 



g(A)-P+^'(g(A)) g(A)-P+^'(g(A)) 



g{A)a p g(A) 1 = g(A)v A>p . 



□ 



2.1 Remark. If a p > 0, d — 1 and f?(A) = 1, then the unique vectors u A and defined 



in (ii) of Proposition 12. II take the forms u A = ~l p with l p :- 



v A 



V 



cl\ + 2a 2 H + pap 



eP xl , and 



a 2 + \- a p 



□ 

2.2 Definition. An INAR(p) process (A n )n>-p+i wii/i coefficients called prim- 

itive if 

(i) a p > ; 

(ii) d = 1, where d is the greatest common divisor of the set {z G {1, . . . ,p} : a» > 0}. 

2.2 Remark. If a p = and there exists i 6 such that a-j > 0, then (A n ) n ^_ p+1 

is an INAR(p') process with coefficients ai,...,a p r with a p > > 0, where p' = maxji G 
{1, . . . , p} : ai > 0}. If a p > 0, but d ^ 2, then the process takes the form 

X k = a d o X k _ d H h a(p/d-i)d ° ^fc-( P /d-i)d + « P ° A fc _ p + e k , k G N, 

and hence the subsequences (Jfd„_.;)„^_p/d-|_i, j = 0, 1, . . . , d — 1, form independent primitive 
INAR(p/d) processes with coefficients ay, «2d, • • • , ot p such that X^ p+d -j = X- p+2d ~j = ■ ■ ■ = 
X_j = 0. Note also that in this case not all of the coefficients ad, «2d, ■ ■ ■ ,a p are necessarily 
positive. Finally, we remark that an INAR(p) process (X n ) n ^_ p+1 is primitive if and only if 
its matrix A defined in (I2.2p is primitive. Indeed, if (A„) n j,_ p+ i is primitive, then part (ii) of 
Proposition [2J] readily yields that A is primitive. Conversely (using the notations of the proof 
of Proposition 12. ip . if A is primitive, then, by the proof of part (i) of Proposition 12. 1\ the 
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digraph D(A) is strongly connected. This yields that a p > 0, since otherwise there would 
be no path from 1 to p. Further, the primitivity of A yields that the index of imprimitivity 
of A equals 1. Using that the cycles of D(A) are 1— > i — > (i — 1) — > 2 — > 1 for all 
i — 1, . . . ,p such that a» > (not considering rotations) and such a cycle has length i, we 
get d = 1. □ 

The next proposition is about the limit behavior of E(X&) as k — > oo. This proposition 
can also be considered as a motivation for the classification of INAR(p) processes, see later on. 

2.2 Proposition. Let (X n ) n ^_ p+ i be an INAR(p) process such that Xq = X_i = • • • = 

X_ p+1 = and E(e^) < oo. Then the following assertions hold: 



(i) // q(A) < I, then 



lim E(X fc ) - 



(ii) // q(A) = 1, then 



k^oc i - Y7 i= i a i 



lim // 1 K(.V ; ! 



fc^oo * " ' Y7i=l ia i ^'(i)' 

where ip is the characteristic polynomial of the matrix A defined in 

(iii) If q(A) > 1, then 



kdT? tv \ d ^e _ dfi e g{A)P 1 



hm g{A)- E(X^) { ^ A)d _ i} kakQ{A) - k {Q{A)d _ l Me{A )) 

for all j — 0, 1, . . . , d — 1, where d is the greatest common divisor of the set [i e 
{l,...,p} : a { > 0}. 

Proof. If «! = ■■■ = a p = 0, then f?(A) = and X& = e k , k £ N, which yields that 
linn^oo E(Xfc) = /i £ , i.e., part (i) is satisfied in the case of a± = ■ • • = a p = 0. If not all of 
the coefficients ati,...,a p are 0, then, by Remark 12.21 (X n ) n j>_ p+ i is an INAR(p') process 
where p' = max {z G {1, . . . ,p} : aj > 0}. Hence in what follows we may and do suppose that 
the original process (X„) n ^_ p+1 is such that a p > 0. 

First we prove the proposition in the case of a p > and d — 1, i.e., in the case of 
(X„) n> _ p+1 is primitive. 

Proof of (i) in the case of a p > and d = 1: In this case we verify that 

oo 

lim E(X fe ) = fj, s ejY i A j e 1 = fi £ ej(l p - A)~ 1 e 1 = ^ . 

fc^oo * — ' 1 — oji — • • • — a„ 

3=0 1 P 

By 02.4)1 . it is enough to prove that if g(A) < 1, then the series Y^LqA 3 is convergent and 
its sum is (J p — A)' 1 . By 02.9)1 . we have 

oo oo oo oo 

]T || < £>W {\\Q(A)-iAi - n A || + ||n A ||) ^^QiAYcAr^ + ^ e (Ay\\n A \\ < oo, 

i=0 3=0 3=0 3=0 
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since g(A) < 1 and ta < 1- One can give another proof for the convergence of E^lo ll^ll- 
Indeed, by Horn and Johnson Corollary 5.6.14], we have g(A) = lim^^ ||A"|| 1//n and 
hence comparison test yields the assertion. Finally, by Lemma 5.6.10 and Corollary 5.6.16 in 
Horn and Johnson [26J, we have Y^jLo-A 3 = (J p — A)^ 1 , and hence, by Cramer's rule, 

eJ(I p -A)- 1 e 1 



det(J p — A) ip(l) 1 — — — a p 
Proof of (ii) in the case of a p > and d — 1: In this case we verify that 



lim k 1 E(X fc ) = [i £ eJn A e x 



He He 



By Q23}, we get 

fe-i fe-i 
E(X fe ) = /i^ ^ A J ei = fi e ej i u A + (A* - n A )) ei 

3=0 j=0 

fe-1 

= k^ £ eJli A e 1 + ^t E e7 ^(^ - n A )e 1; fc e N. 

3=0 

By (ESD, we have 

oo oo 

^ || A j - n A || < ^c A r A < oo, 

i=o 

where denotes the pxp nullmatrix. This implies lim^^ k~ x E(X^) = fi £ ejTI A e±. By 
Proposition 12 .![ in the case of a p > and d — 1 (f?(A) is not necessarily 1) we have 

(2.10) e a n A ei = ei u A . A ex = u A>1 v A>1 = g{Ayp+1 ^ g{A)) = ^^y- 

By f )2.7p . we have ax + h a p = 1, and hence 

1 1 



which yields that 



e l Il A ex 



y?'(l) p- (p- l)«i - (p - 2)a 2 2a p _ 2 - « P -x 

1 



EL i 1 *«i + P (i - ELi 1 o<) ' 
which yields part (ii) in the case of a p > and d — 1. 
Proof of (iii) in the case of « p > and d = 1: In this case we verify that 

lim g(A)- k E(X k ) = J£ eJU Aei - 



" WA)-iy(^))' 



ii 



By (EQ|, we get for all k G N, 



fe-i 



fc-i 



^(A)" fc E(X fc ) = e (A)-V e e7 J] A* ei = g(A)- k fi £ eJ ^ (£?(A)->n A + (A? - e (A)^n A ))ei 



j=0 



fe-i 



3=0 
fe-1 



^ej ^ e (i4)^ fc n A Ci + ^eJe(A)- fc ^(^" - f?(A)''n A ) ei . 

j'=o j=0 



Since £?(A) 1 < 1, we have 



fe— 1 k 
3=0 1=1 



1 



f?(A) - 1 



as — ?• oo. 



Further, by (123) . for all k G N, 



fe-i 



e(A)-*53(A^- ff (Ayn x ; 

i=o 

If £>(A)r A 7^ 1, then 



fe-i fe-i 

< 53ff(A)-^|| ff (A)-^ - n A || ^ c A ^(A)- fc +^- 

i=o i=o 



fc-i 



j=0 



1 - g(A)r A 

since £>(A) > 1 and r A < 1. If ^(A)r A = 1, then 

fc 



< / ^ ^ -> as k ->■ oo, 



fc-i 



^(A)- fc ^(A^-^(Ayn A ; 

i=o 



e(A)' 



— >• as — > oo. 



Using also (12.61) and (12.101) . this concludes (iii) in the case of a p > and d — 1. 

Now we turn to give a proof in the case of a p > and d ^ 2. In this case, by Proposition 
12. 1[ A is irreducible, £?(A) > and, by Remark [2T2| the subsequences (Xd,n-j)n^-p/d+u 3 — 
0, 1, ... , d—1, form independent primitive INAR(p/d) processes with coefficients a d , aid, ■ ■ ■ ,a v 
such that X_ p+ d-j — X_ p+ 2d-j = ■ • 



X_j = 0. Let us introduce the matrix 



1 

1 













1 



l(p/d)x(p/d) 



and its characteristic polynomial 

£(A) := det(A/ p/d — A) = A p / d - a^- 1 - a 2d X^ d ~ 



- a p _ d \ - a p , A G C. 



Since the greatest common divisor of the set {i G {1, . . . ,p/d} : aid > 0} is 1, by Proposition 
I2.1[ we have A is primitive. We check that g(A) d = g(A). Since <f(X) = <£>(A d ), A G C, 
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we get g(A) d ^ f?(A). By Proposition 12.11 g(A) > and f?(A) is an eigenvalue of A. 
Hence g(A) l l d is an eigenvalue of A, which implies that f?(A) ^ g(A) l l d or equivalently 
g(A) d > g(A). 

If g(A) < 1, then g(A) < 1 and using that part (i) has already been proved for primitive 
matrices (i.e., in the case of a p > and d = 1) we have for all j = 0, 1, . . . , d — 1, 

lim E(X nd / /? 



1 - a d - a 2d ct v 1 - Yh=i a i 

This yields that lim^^ E(X n ) exists with the given limit in (i). 

If g{A) = 1, then g(A) = 1 and using that part (ii) has already been proved for primitive 
matrices we have for all j — 0, 1, . . . , d — 1, 

E(X dn ^ j ) _ /j £ d\i e dfx £ 



lim 



n a d + 2a 2d H h ~ d a v da d + 2da 2 d H h pa P XlLi ^ 



This yields that lim^oo k 1 E(Xfc) exists with given limit in (ii). 

If g(A) > 1, then £>(A) > 1 and using that part (hi) has already been proved for primitive 
matrices we have for all j — 0, 1, . . . , d — 1, 

j. E(^rad-j) _ /4 _ ^£ 

™ e (A)» (^(A) - 1) Yl=i kotkdg(A)- k ~ (g(A) d - 1) Ylk=i kda kd g(A)- kd 
_ dfe dp, e g(A) p ~ 1 

~ (g(A) d - 1) EU^( A )- £ ~ (e( A ) d - 

where the last equality follows by (I2.6j) . Since 

E(X fc(J _j) E(Xfcd_j) 
hm . ...j = hm ~ , j = 0, 1, a - 1, 

k^oo g{A) kd k^oo g[A) k 



we have (hi). □ 



Based on the asymptotic behavior of E(X^) as k — > oo described in Proposition 
we distinguish three cases. The case g(A) < 1 is called stable or asymptotically stationary, 
whereas the cases g{A) = 1 and g(A) > 1 are called unstable and explosive, respectively. 
Note also that, if a p > 0, then, by ( 12. 71) of Proposition 12.11 g(A) < 1, f?(A) = 1 and 

g(A) > 1 are equivalent with a± + h a p < 1, oti + • ■ • + a p = 1 and a\ + • — h a p > 1, 

respectively. 



3 Convergence of unstable INAR(p) processes 

A function / : R + — )■ R is called cddldg if it is right continuous with left limits. Let D(R + , R) 
and C(R + ,R) denote the space of all real- valued cadlag and continuous functions on R + , 
respectively. Let T)^ denote the Borel ex-field in D(R + ,R) for the metric defined in (16.4) 
in Billingsley [5] (with this metric D(R+,R) is a complete and separable metric space). For 
stochastic processes (!Vt)teiR+ and (!Vf)teiR +) n e N, with cadlag paths we write y n — > y 
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if the distribution of y n on the space (D(R + , R), T> 00 ) converges weakly to the distribution 
of y on the space (D(R + , R), V^) as n — > oo. 

For each neN, consider the random step processes 

X? := rrtf^j, t G R + , nGN, 

where [^J denotes the integer part of a real number x G R. The positive part of x G R will 
be denoted by x + . 

3.1 Theorem. Let (Xk)^_ p+ i be a primitive INAR(p) process with coefficients oti, . . . , ot p G 
[0, 1] swc/i i/iai ai + • • • + a p — 1 (hence it is unstable). Suppose that X = X-\ = ■ ■ ■ = 
X_ p+ i = and E(ef) < oo. Then 

(3.1) X n ^ X as n^r oo, 

where (X t )t^R + is the unique strong solution of the stochastic differential equation (SDE) 

(3.2) dXt = J—(ji e dt + y/ alX? dW t ) , t G R+, 

with initial value X = 0, where 

cp'(l) = a\ + 2« 2 H Ypct v > 0, a 2 a := «i(l - «i) H h a p (l - ot p ), 

and (Wt)teR + is a standard Wiener process. (Here <p is the characteristic polynomial of the 
matrix A defined in (12.21) . ) 

3.1 Remark. Note that under the conditions Theorem I3.1[ if p ^ 2, then er^ > 0, and if 
p — 1, then cr^ = 0. Indeed, if p ^ 2, then a p < 1, since otherwise «! = ••• = = 
and hence the greatest common divisor of {i G {1, . . . ,p} : o.i > 0} = {p} would be p, which 
is a contradiction. Since, by our assumption a p > 0, we get <j\ ^ a p (l — a p ) > 0. If p = 1, 
then a p = ai = 1, and hence a\ = ai(l — «i) = 0. 

Remark also that in the case of p = 1 we have a± = 1 and hence X n = 5^=1 71 ^ ^ 
(p'(l) = 1, o"^ = and then the limit process in Theorem 13 .11 is deterministic, namely X t = /j, e t, 
t G R + . To describe the asymptotic behavior of an unstable INAR(l) process one has to go 
one step further and one has to investigate the fluctuation limit. By Donsker's theorem (see, 
e.g., Billingsley [5j Theorem 8.2]), we have ^/n(X n — E(X n )) a £ W as n — > oo, where W 
is a standard Wiener process. For completeness, we remark that Ispany, Pap and Zuijlen [2"7"| 
Proposition 4.1] describes the fluctuation limit behavior of nearly unstable INAR(l) processes. 
□ 

3.2 Remark. The SDE (13 .2p has a unique strong solution (Xf) t ^o for all initial values 
Xq = x G R. Indeed, since \y/x — ^/y\ ^ \J\x — y\, x,y ^ 0, the coefficient functions 
R 3 x H- /i e /<//(l) and R 3 x h> ^ a^x + / \p' (1) satisfy conditions of part (ii) of Theorem 
3.5 in Chapter IX in Revuz and Yor |52j or the conditions of Proposition 5.2.13 in Karatzas 
and Shreve [37] . Further, by the comparison theorem (see, e.g., Revuz and Yor [52} Theorem 
3.7, Chapter IX]), if the initial value Xq = x is nonnegative, then X* is nonnegative for 
all t G R_|_ with probability one. Hence Xff may be replaced by X t under the square root 
in (13. 2p . The unique strong solution of the SDE (13.21) is known as a squared Bessel process, a 
squared-root process or a Cox-Ingersoll-Ross (CIR) process. □ 
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3.3 Remark. If the matrix A is not primitive but unstable, then we can suppose that a v > 0, 
since otherwise it is an unstable INAR(p') process with p' := max {i G {1, . . . ,p} '■ a» > Oj 
(note that there exists an i G {1, . . . ,p} such that on, > because of the unstability of A). If 
a p > and d ^ 2, then, by Remark [2T2| the subsequences (Xrf n _j) n ^_ p /rf + i, j — 0,1, ... , d—1, 
form independent primitive INAR(p/d) processes with coefficients a^, a 2( 2, . . . , a p such that 
A-p+d-j = X_ p+2 d-j = ■ ■ ■ = X^j = 0. Hence one can use Theorem 13 .11 for these subsequences. 
With the notations 



n 



t G R + , + J =0,1, 

d 



d-1. 



by Theorem 13.11 X n,:1 
of the SDE 



ax; 



X^> as n _ ). 00) where (A^) te R + is the unique strong solution 
1 



a d + 2a 2 d H h ^oi. 



;(?) 



t G 



with initial value Xq = and (W t )t e K + , j = 0, 1, . . . , d — 1, are independent standard 
Wiener processes. We note that if a p > and d ^ 2, then Af n does not converge in general 
as n — > oo. By giving a counterexample, we show that even the 2-dimensional distributions 
do not converge in general. Let p := 4, ot\ = «3 := 0, a 2 = «4 := 1/2. Then d = 2 and 
using that X n ^ 



— > X^ as n — > 00, j = 0, 1, we have 



(3.3) [X^'°,X^° 
as n — > 00, and 

(3.4) [A^ 1 ,*; 



— A2„, — X^n 

n n 



(o) v (o)i 



converges in distribution to [X{ , X 2 



— A2„_l, — X4„_i 



(i) vWl 



converges in distribution to \X{ , X : 



as n — > 00, where (A^) t6 R, is the unique strong solution of the SDE 



d^' } = ~ dt + J^o*f y dn^) , t g 



with initial value <Y = 0, j = 0, 1. However, we show that 



(3.5) [X?,X< 



Indeed, we have 



— X n , —X 2n 
n n 



\xl n ,xl n \ 



does not converge in distribution as n — > 00. 



x x 1 X ' 

2n 2n 



and hence, by (13. 3ft . 

[Af 2n , Af 2 2n ] converges in distribution to 



1 -v n )0 v n )0 

2 1 2 2 



1 v(0) 1 «(o) 
2^i '2^ 2 



as n — >• 00. 
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Further, using that 



[x? n - 1 ,x* n - 1 ] 



2n - 1 



-X 



1 



2n-lj 



2n - 1 



-X 



2(2n-l) 



n 



n,l v 2ra-l,0 



2n - 1 



^1 , <*i 



and that the subsequences (X 2n _i) ne rc and (X 2 (2 n -i))neN are independent, by (13. 3p and (13 .4p . 
we get 



[X^ n 1 ,A' 2 2n x ] converges in distribution to 
Since the random variables 

1 ,„s 1 

and 



1 v(0) ^(O) 

2^i 



■v(l) v(0) 



v(0) 

2^i 



as n — y oo. 



do not have the same distributions (the coordinates of the first one are dependent, however the 
coordinates of the second one are independent), we get (13.51) . □ 



For proving Theorem 13 .1[ let us introduce the sequence 

(3.6) M k :=X k -E(X k \F k _ 1 )=X k -a 1 X k „ 1 a p X k ^ p - /i £ , fcGN, 

of martingale differences with respect to the filtration (J-fc)fcez+, and the random step processes 



[nt\ 



M?:=n- iy %2M k , t 



n e N. 



k=l 



First we will verify convergence 

(3.7) M n A M asii4 oo, 

where (M.t)teR+ is the unique strong solution of the SDE 



(3.8) 



dM t 



o'(l) 



with initial value .Mo = 0. The proof of (13.71) can be found in Section HI 

3.4 Remark. If (Xf) te m. + is a strong solution of (13. 2p with initial value X£ — x e M, then, 
by Ito's formula, := ip'(l)X* — /i e t, t G R+, is a strong solution of (13. 8p with initial value 
JHq = <f'(l)x. On the other hand, if {Ai\)tm+ is a strong solution of (13.81) with initial value 
Ml = y E R, then, again by Ito's formula, 



(3.9) 



t e 



is a strong solution of ( 13. 2ft with initial value Xq = -^rmV- Hence, by Remark 13. 2\ the SDE 
( 13. 8 p has a unique strong solution (Ai y )t>o for all initial values A^q = y G M. Further, if the 
initial value M. v = y is nonnegative, then Ai v t is nonnegative for all t e M + with probability 
one. Hence (.M t + fi £ t) + may be replaced by (Ait + under the square root in (I3.8p . □ 
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Moreover, from ( 13.61) we obtain the recursion 

X k = a x X k _i H h a p X k _ p + M k + /i e , k e N, 

which can be written in the form X k = AX k ^\ + (M k + // e )ei, fceN. Consequently, 



implying 
(3.10) 



X k = J2( M J + »s)A k - j ei , fceN, 
i=i 



X fc = e[X fc = ^(M,- + fi £ )ejA k -i ei , keN. 

3=1 



In Section HJ we show that the statement (13.11) will follow from (13. 7p and (13.101) using a version 
of the continuous mapping theorem (see Appendix). 



4 Proof of Theorem 13.1 



For the proof we will use Corollary 15. 1[ Theorem 15.11 and Lemma 15.21 which can be found in 
Appendix. 

______ £ 

First we prove (IXTj) . i.e., M n ^ M as n — > oo. We will apply Theorem 15 .11 for U = Ai, 
U* = n~ l M k) n, k en, and for {J r k ) k&+ = {J r k)k&+, n G N. By Remark E3 the SDE 
(13.81) has a unique strong solution for all initial values Ai^ = x, Now we show that 

conditions (i) and (ii) of Theorem 15.11 hold. We have to check that for each T > 0, 



(4.1) 
(4.2) 



sup 

te[o,T] 

\nT\ 



™ ~ V (1) Jo 



J_ E(M fc 2 l { , A4 | >ne} | .F fc _i) -A for all 9 > 



as n — >■ oo, where 



fc=i 
p 



means convergence in probability. 



By ( 13. 6\\ and using also that at\ 



a p = 1, we get 



|rasj 



M n s + fl e S = U [ X k aiXk -i _ I -r /' . s 

fe=l \ i=l 



p— i L ns J — * 



// | A fc -^ai X fe ] + 

,fc=LnsJ— p+1 »=1 fc=[nsj— p+1 



7? 



/'c- 



p p 



n 



He- 



ll 



Thus (M™ + /i e s) + = + fi £ s, and using that 

' ns - [ns\ t 2 111 L? ^ 1 / [nt\ \ \ [nt\ + (nt - [nt\f 



n 2 n \ n ' V n / / 2n 2 

y fe=i 

we get 

/"* 1 LntJ_1 P P - I nil P P 

/ (m: + fi £ s) + ds = - EE a ^+ i+ — d^EE a ^j-i+i 

JU e=o j=i i=j j=i i=j 



[nt\ + (nt - [nt])'" 
2n 2 



H TTt V*- 



Hence, using that <p'(l) — o.\ + 2a2 + ■ • • + pa p , we have 

[ntj-l 1 p \nt\-l 



r t , m [ntj-l p [nt\-l 

/ (m:+w)+6s = *£>- y: Xi-^T,* E * 

I rt lb lb 



(-=0 i=2 j = l_ nt \-i+l 

nt - [nt\ v [nt\ + (nt - \nt\ ) 2 

+ ~2 1^ a i X lnt}- j+ l + ^ 



Using (15.41) . we obtain 



fc=l fc=l \i=l / 

I P f [n*H+l \ 

^E «*(!-«*) E ^i-O + 



[nt] 2 



, > . , - ■ 9 £ 

rz^ 11 — ' \ 14 — ' / n 



^E*--^E ha-" 1 ) E *i +^.- 



Hence, for all neN, the randomness of the difference in (14.11) is via a linear combination of 
the random variables X\ nt \_j, j = 1, . . . , p. Then, in order to show (I4.ip . it suffices to prove 

1 P 

(4.3) sup — g X\nt] — > as n — > oo. 

te[o,T] n 



By f l3~T0D and f O]) . 

\nt\ . \nt\ . 

i=i ^ i=i ' 



Consequently, in order to prove (14. 3p . it suffices to show 

[nT\ 

n l 



— 2J — > as n -)■ 00. 
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n 2 



In fact, one can show that n 2 E(|Mj|) — > 0. Indeed, Corollary 15 . 1 1 yields that 



^E(|M,-|)< — ^VT^— L^Jv / L^I^0 asrwoc, 



z — ' n 

i=i 3=1 

with some constant X G R+. Thus we obtain (I4.ip . 

To prove (I4.2p . consider the decomposition M k = N k + (e k — /i e ), where, by (15. 7p . 

Nk ■= ^2 (%k,i,e - E(^fc,i,^)) H h 2^ (Cfc,p^ _ E(£fc )P ^)). 

£=1 £=1 

Clearly, 

2 /" 2 / \2\ 

^fc ^ 2 (^fc + - A*e) J and l{|M fe |>n0} < l{|7V fc |>n0/2} + l{|e fe -/te|>7i0/2}, 

and hence (14 .2p will be proved once we show 

\nT\ 



(4.4) E(N 2 l {lNk \ >n9} | -A for all 6 > 0, 

fe=i 

(4.5) — ^ E(N 2 l {]£k ^ l>ne} | 7U) A for all 6 > 0, 



fc=i 

LnTj 



(4-6) J_^E(( £fc - /U£ )) 2 |J- fe _ 1 



Ao. 



k=l 



First we prove ( 14.41) . Using that the random variables {£,k,i,j : j G N, i G {1, . . . ,p}} are 
independent of the cx-algebra J^k-i for all k G N, we get 

E(iV fc 2 l { | iV , >n0} | J^k-l) — Fk(Xk-l, • • • , Xfc_ p ), 

where F k : — )■ R is given by 

F fc (zi, . . . , 2; p ) := E((Sfc(2i, . . . , z p ) l{|5 fe ( 2l ,..., 2p )|>ne})), 
with S'fc^i, . . . , z p ) := Y7i=i I2eLi(£,k,i,e - E (&.,*,£)). Consider the decomposition 

F k (z 1: . . . , z p ) = A k (zi, . . . , Zp) + B k (zi, . . . ,Zp), 

where 

A k (z 1 , ...,Z p ) := E((£ k ,iJ - Fj(£k,il)) 2 ~ft-{\S k (z 1 ,...,z p )\>nd}), 

i=i e=i 

B k (zi, ...,Z p ) := E((^fe,i,£ - E{£k,i,t))(£k,j,t> - ^{^k,j,e'))~i-{\S k (zi,...,z p )\>ne}), 

where the sum Yl' is taken for i, j — 1, . . . ,p, £ — 1, . . . , Zi, £' = 1, . . . , Zj with (i, £) ^ (j, £'). 
Consider the decompositions 



S k {Zi, ...,Z p ) — (^k,i,e — E(^k,i,i}) + Sfc,l( Zl i ' ' ' ' Z p)> ^ — 1, . . . ,p, £ — 1, . . . , Zi 
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where 

Sl,e( z ii ■ ■ ■ i z p) '■= ^ faj'*' ~ E (£w))> 
where the sum Yl" is taken for j = 1, . . . ,p and £' = 1, . . . , Zj with (J, £') ^ (i, 
Using that 

l{|5 fe (z 1 ,...,z p )|>n0} < l{|a,<,i-E(e k ,i^)|>nfl/2} + l{|S*^( 2l ,..., Zp )| >7l e/2}' 

we have 



A k (zi, . . .,z p ) ^ A^. 1 (z^ . . . , z p ) + A^\zi, . . . , 2 p ), 



where 



P H 



A% (zi, ...,z p ):=^2^2 E ((£fc,M - E(^ >i/ )) 2 l { | 5fci ^_E( ?fc ^,)|>„e/2}), 

i=l 1=1 

P Zi 

Yl Yl E (fe^ - E ^ k ^ 2l {\Si i (z 1 ,...,z p )\>ne/2})- 



A k \zi, . . . ,Zp) :— 

i=l £=1 

In order to prove (I4.4p . it is enough to show that 

\nT\ _„T\ 



^2 Y A k\ X k~l, X k-p) -A 0, ^2 S ^1 2) (^-1' • • • ' X k 



k=l 



k=l 



1 L«TJ 



fc=i 



(4.7) 

as n — > oo. We have 

p 

^(z!,...,^) = 5^^E((^i ii) i-E(^ ul )) 2 l { | &ilil _ EKl>iil) | >nfl/2 }), fceN, 



i=l 



where 

E ((&M,1 _ E (6,i,l)) 2 l{|6,i,i-E(ei,i, 1 )|>n0/2}) -> 0, 

as n — 7- oo for all i e {1, . . . ,p} by the dominated convergence theorem. Thus, by Corollary 
I5.lt we get with some constant K 6 M + , 

LnTj 



rr 



fe=i 



LnTj p 



p 

i=i 



72 /^ E (^fc-i)E((Cl,j,l - E(^l i i i l)) 2 l{| ?1 ^ 1 _ E (g lilil )|>ne/2}) 



fc=l i=l 



LnTj 



E((^,i,l-E(^i )i) i)) 2 l{| &fiil _E(f 1A1 )|>ntf/2})^ 5^ ^ 



fe=i+l 



LnTj ([nTj + 1) 
2n 2 



2jE((^i jij i - E(^i i i i i)) 2 l { | 6 ^ 1 _ E{6 ^ l) | >ne/2 }) ->> 0, 
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which yields n~ 2 £)H A? (X k _ u . . . , X fc _ p ) A 0. 

Independence of — E(£ fc ,i,£) and S l kt (zi, . . . , z p ) implies 



p z% 

I (2), 



4 „ , 4 w , . 4 

i=i 



4 p p 



4 2) (z 1; . . . , z P ) = ee e ^ - E (^)) a ) pd^i. • • • > *oi > ^/ 2 )- 

Here E((£ fc)i ^ — E(£ fci ^)) 2 ) = a^(l — i = 1, . . . , p, and, by Markov's inequality, 
P(|^i, - - - , z p )\ > n9/2) <: -L E(Si 4 ( Zl , z p f) 

tv v ■ Vax(^(zi, ...,zp)) = ^ E" a ^ " a ^ ^ ^ E z o a oi l ~ 

Thus we get 

A { k 2) (z u ...,z p )^ E E z i z i a ^ 1 ~ a i) a A l ~ a j 

i=l j=l 

Hence, by Cauchy-Schwarz's inequality and Corollary 15. 1[ we get with some constant K e R + , 

\nT\ \nT\ p p 

- E E(4 2) (^-i: ■ • • , **-p)) < E E E E^^X^Oa^l - - a,-) 

fc=i fc=i i=i j=i 

^E fc '£<*(!-«*) -+°> 
fc=i \i=i / 

which implies n~ 2 Y}k=i A [ 2 \ x k-i, ■ ■ ■ , X k _ p ) — >■ 0. 
By Cauchy-Schwarz's inequality, 



l^fcOi, •••,%) | < yj B^izu ...,z p ) E(l { | Sfc(zii ... i2p) | >ne} ), 

where 

Bjp(z u ...,Zp) := E f(E _ E(^))(Cfcj,£' - E(Cfe,^'))) Y zi,...,^, G Z+. 

Using the independence of 6fe,M _ E (6m,*) and £ fcj -^ - E(£ fc> ^>) for (i,£) ^ (j,£'), we get 

Bj^Czi, . . . , Zp) = E Oi(l - ai)^'! 1 - a j) 
p 

= Y Z i( Z i ~ l ) a2 i( l ~ a if + E Z i Z i a i( l ~ ai)aj(l - 0£j) 

1=1 

^Ki(zi + --- + 2; p ) 2 , 
with some constant K\ e R + . Further, by Markov's inequality, 

1 P if 

E(i { | 5fc( , 1 ,..., 2p) | >ne} ) < e - a i) ^ ^ (zi + ' ' ' + Zp] ' 

i=i 
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with some constant K<i G R+. Hence 

K 

\B k (z u ...,z p )\ ^ — {zx H h z p f /2 , zx,...,z p e Z+, 

n 

with some constant K G K+. Using that 

(*! + ■■• + z p ) 3 / 2 < c p (^ /2 + ■ ■ • + ^ /2 ), Zl,..,^£Z + , 

with some constant c p G R+, we get, in order to show (14.71) . it suffices to prove 
n~ 3 ES (*22 + • • • + Xl%) A 0. In fact^ n"» EH ( E(X^) + ■ ■ • + E(X 3 / 2 p) ) 
since Corollary 15.11 implies E(X 3 ^ 2 ) ^ (E(X|)) 3//4 = 0(£ 3 ^ 2 ). Thus we finished the proof of 

(33). 

Now we turn to prove (14. 5p . Using that for all k G N the random variables {£/c,i,i> £fc : J G 
N, i G {1, . . . , p}} are independent of the cx-algebra J-fc_i, we get E(N^t^ £k _^\ >ne y \ J-u-i) = 
G fc (Jf fc _i, . . . , X fe _ p ), where G fc : %\ ->• R is given by 

Gk(zx, . . . ,Zp) := E(5'fc(2;i, . . . , Zp) 2 l{| efc _ iUe |> ri e}), 21, . . . , 2^ G Z + . 

Using again the independence of {^jfe,ij, £fc : j G N, z G {1, . . . ,p}}, 

G k (z u ...,z p ) = P{\e k - f ji e \> n9) £ J] E((£ fe ,^ - E(£ fc , M )) 2 ), 

8=1 1=1 

where by Markov's inequality, P(|ejfe — )U e | > n9) < rr 2 Q~ 2 E((e k — fi £ ) 2 ) = n~ 2 6~ 2 a 2 , 
and E((^ fc j£ — E^,^)) 2 ) = — Hence, in order to show (14.51) . it suffices to prove 

n" 4 EH ^fc 0. In fact, by Corollary EU n~ A E ( X k) -> 0. 

Now we turn to prove (I4.6p . By independence of and J-^-i, 

1 [nTJ 1 LnTJ I Tl 

- £ E((e k - ^) 2 1 7- fe _0 = -J2 E ((^ " ^) 2 ) = ^T°l ^ °> 
n z * — ' n z *■ — 4 n z 

k=l k=l 

thus we obtain ( 14. 6p . Hence we get (14. 2p . and we conclude, by Theorem 15. 1[ convergence 

Mn^M. 

Now we start to prove (13. ip . By (I3.10p . X n = ^/ n (A4 n ), where the mapping \l/„ : 
D(R + ,R) -»■ D(R + ,R) is given by 

^(/)W:=|(/(g-/(^) + ^)erA«^ 

for / G D(R+,R), t G R+, n G N. Further, # = where, by (133]) . the mapping 

* : D(R + ,R) -> D(R + ,R) is given by 

*(/)(*) :=-^t(/(*)+^). /gD(R + ,R), t G R+. 

We check that the mappings n£N, and ^ are measurable. Continuity of ^ follows from 
the characterization of convergence in D(R + ,R), see, e.g., Ethier and Kurtz [181 Proposition 
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3.5.3], thus we obtain measurability of tf. Indeed, if f n G D(R+,R), n G N, / G D(R+,R) 
and the sequence (f n )w=N converges in D(R + ,R) to /, then for all T > there exist 
continuous, increasing mappings A n , n6N, from [0, oo) onto [0, oo) such that 



lim sup | \ n (t) — t\ 

n ^°°te[o,T} 



Since for all t G R + 

|¥(/„)(A„(f)) -*(/)(*)! 







and lim sup |/„(A n (t)) - /(t)| = 0. 

n ^°°te[o,T] 



1 



^(1) 
1 



C/n(A„(t))+/i e A n (t)) 



^(1) 



(/(*) + A*.*) 



/4 



-T7TTl/n(A n (t))-/(t)| + -^-|A n (t) 



we have for all T > 0, 



lim sup mf n )(\ n {t))-^(f)(t)\ = 0. 

n ^°°tG[0,T] 



In order to prove measurability of \& n , first we localize it. For each N e N, consider the 
stopped mapping V% : D(R + ,R) -»• D(R + ,R) given by **(/)(*):=*„(/)(* A JV) for 
/GD(R + ,R), £ G R + , n,JVGN. Obviously, *^(/) -»• * n (/) in D(R + ,R) as iV^oo for 
all / G D(R + ,R), since for all T > and N^T we have := *«(/)(*), * G [0,T], 

and hence sup te r 0T i (/)(£) — ^n(/)(*)| - >• as AT — > oo. Consequently, it suffices to show 
measurability of for all n, N G N. We can write = ^/^' 2 o Sl/^' 1 , where the mappings 



^ : D( 



M->- 



and ^' 2 



;m-/V+l 



->• D( 



|n(tA7V)J 



are defined by 



^' 2 (X ,X!, 



,^)(t):= £ (xj-x^ + ^ejAM-tei 
i=i n 

for / G D(R + , R), t G R+, x = (x , x niV ) G R nAf+1 , n.iVeff. Measurability of ^ 

follows from Ethier and Kurtz [T8| Proposition 3.7.1]. Next we show continuity of ^n' 2 by 
checking sup tg [ 0T ] \^f^ ,2 (x k )(t) — ^/^ ,2 (x)(t)\ — > as k — )■ oo for all T>0 whenever :r fc — )■ x 
in R nAr+1 . This convergence follows from the estimates 

|n(TA7V)J 

sup |v^V)W - *n' 2 (*)WI < E (14 - x il + - c3"AL«*J-i Cl 
te[o,T] j=1 

since ejA^^e-y < ||A Ln * J_J '|| ^ C A . We obtain measurability of both and ^- 2 , 

hence we conclude measurability of SE^. The aim of the following discussion is to show that 
there exists C C C*,(*„)„ 6N with C G and P(.M G C) = 1, where C*,(*„) nef) is defined 
in Appendix. We check that C := {/ G C(R + ,R) : /(0) = 0} satisfies the above mentioned 
conditions. First note that C = C(R + , R) n 7i \0), where vr : D(R + ,R) -> R, tt (/) := /(0), 
/ G D(R + ,R). Using that C(R + ,R) is a measurable subset of D(R + ,R) (see, e.g., Ethier 
and Kurtz [TBI Problem 3.11.25]) and that 7r is measurable (see, e.g., Ethier and Kurtz [T8| 
Proposition 3.7.1]), we have C G V^. Fix a function / G C(R + ,R) and a sequence (/ n )neN 
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in D(R + ,R) with f n — ^> /, where is defined in Appendix. By the definition of we 



get *&(f) G C(E + ,M.). Further, we can write 

' 3 -1\ , fie 



\nt\ 
5=1 ^ 



/n 



n 



e 1 n A ei 



5=1 



n 



fn 



1 ' ^+^)eJ(A^-n A )e 1 , te 



Using (I2.10p and the assumption g(A) — a± + h a p — 1, we get ejYlA e i = and 



5=1 



n / n 



- =/« 



/ B (o) + M^. 

n 



Thus we have 



l*n(/»)W-*(/)WI<-^; 



[nt\ 
5=1 



f n 



fn 



[nt\ 



n 



-fit) 

fn 



_|_ ^ + l/« 



n^(l) ^'(1) 



n 



Here for all T > and t G [0,T], 



fn 



n 



fit) 



[nt\ 



n 



-f 



[nt\ 



n 



f 



[nt\ 



n 



-fit) 



^ L0T (f 1 n- 1 )+ sup |/ n (*)-/(*)|, 



te[o,T] 



where wr(/, •) is the modulus of continuity of / on [0, T], and we have wt(/, k" 1 ) —> 
since / is continuous (see, e.g., Jacod and Shiryaev [301 Chapter VI, 1.6]). In a similar way, 
for all j = 1, . . . , \nt\ , 



fn - - /, 



J-l 



/? 



^ Wr (/,^ 1 ) + 2 sup \f n (t) - f(t)\. 
te[o,T] 



By (EID, since g(A) = 1, 



[ntj [n*J 

- n A || < Y, c ^' j < 7^ 

5=1 ' 



5=1 



Further, 



|/»(0)| < |/„(0) - /(0)| + 1/(0)1 ^ sup \f n (t) - f(t)\ + |/(0)|. 

te[o,T] 



Thus we conclude C C C#,(#„)„ eN . Since Ai = and, by the definition of a strong solution 
(see, e.g., Jacod and Shiryaev [2DJ Definition 2.24, Chapter III]), M. has almost sure continuous 
sample paths, we have P(A4 G C) = 1. Consequently, by Lemma 15.21 we obtain X n = 
^niM n ) -A V(M) = X as n -»■ 00. □ 
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5 Appendix 



In the proof of Theorem 13.11 we will extensively use the following facts about the first and 
second order moments of the sequences (X k ) k£ z + and (M k ) k& z + - 

5.1 Lemma. Let (X k ) k ^^ p+ i be an INAR(p) process defined by ( 12.1 j) such that X = X_i = 
• ■ ■ = X_ p+ i = and E(ej) < oo. Then, for all fceN, 

fc-i 

(5.1) E(X k ) =/i £ J2 e i Aie ^ 

1=0 

k— 1 V k—i—1 j 

(5.2) Var(X fc ) =a e 2 ^(e^e 1 ) 2 + ^^a l (l-a l ) £ J^A^^eO'Ce^eO- 

£=0 i=l j=0 «=0 

Moreover, 

(5.3) E(M fc | JF^) = /orfcGN, 



(5.4) E(M fe M £ | J- max{fc , n . 



ai(l - «i)X fe _! H h a p (l - a p )X fe _ p + a 2 if k = £, 

if k^£. 



Further, 

(5.5) E(M fc ) =0 /or fceN, 



(5.6) E(M fc M, 



ai(l - ati) E(X k ^x) H h a p (l - a p ) E(X fc _ p ) + a 2 if k — 

if k ^ 



Proof. We have already proved (I5.ip . see (12.41) . The equality M k = X k — E(X k | J-^-i) clearly 
implies Q and (1531) . By ((2HJ) and (TO . 

(5.7) M fc = ]T (£ Mj - - E(£ Mii )) + • • • + X) few - E(£ fcjPJ )) + (e* - E(e fc )) • 
i=i i=i 

For all k G N, the random variables {Cfc,ij — E(^ i jj), e k — E(e k ) : j 6 N, i 6 are 
independent of each other, independent of J- k -i, and have zero mean, thus in the case k = I 
we conclude (JEU) and hence flSTBP . If fc < £, then E(M k M e | 7^_i) = M fc E(M/ | = by 

(15.31) . and thus we obtain (15.41) and (15. 6p in the case of k ^ t. 

By (13.101) and (15. ip . we conclude 

k 

X k - E(X fc ) = X M je jA k - j ei , fceN. 

3=1 
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Now, by flEED, (EUD, 

k k k 

Var(X fc ) = ^^V(M j Mz)e T 1 A k - ] e l e T l A k - t e 1 = ^ E(M 2 )(eJ A k - j e x ) 2 
j=i 1=1 j=i 

= E ft"^ 1 - «*) E (*i-«) + ^ 2 ) (cr^ci) 8 
j=i \i=i j 

k p k 

= a 2 J2(ejA k -J ei f + J2 «*(! - «i) E E(X J - l )(e[A^e 1 ) 2 , 
J'=l »=i i=i 

and hence, using also that E(X ) = E(X_i) = • ■ ■ = E(X_ p+ i) = 0, we get 

k— 1 p & 

Var(X fc ) = ( x £ 2 ^(e^e 1 ) 2 + ^a l (l-a i K £ ]T (e^e^^A^eO 2 

£=0 i=l 3=i+l £=0 

= a 2 £>7 A'd) 2 + a*. E - «*) E E( e ^ A ' e i)(^ Afe ^ le i) 2 ' 

£=0 i=l i=0 £=0 

which yields (JOJ- D 

5.1 Corollary. Lei (X k ) k ^_ p+ i &e a primitive INAR(p) process defined by (I2.ip suc/i £/iat 
ai H — • + a p — 1 fz.e. unstable), X = X_i = ■ ■ ■ = X_ p+1 = and E{e\) < oo. Then 

E(X k ) = 0(k), E(X 2 ) = 0(k 2 ), E(\M k \) = Oik 1 / 2 ). 

Proof. By (loTTj) . 

fc-l 

E(X fc )^/i £ ^||A £ |KC A /i e fc, 

where 



(5.8) C A '■— sup < oo. 

Here is finite since, by (12.91) . Ca ^ ca + 1 1 || . Hence we obtain E(Xfc) = 0(k). We 
remark that E(X^) = O(k) is in fact an immediate consequence of part (ii) of Proposition 12.21 

We have, by Lyapunov's inequality, 

i 1/2 

txiyi — uci) myv\ k _i) ~r u 2 

i=l 

/ p \ 1/2 

hence we obtain E(|M fc |) = 0(k 1/2 ) from E(X fc ) = 0(k). 
Thus we get 

E(X 2 ) = Var(X fe ) + (E(X fe )) 2 = 0(k 2 ). 



E(\M k \) <: JE(M 2 ) = 53 0,(1 - «,) E(X fc _ 4 ) + a 2 
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Indeed, by fl£2]) and (EH 



e=o i=i j=o 1=0 



where cx^, is defined in Theorem 13.11 Hence we obtain E(X 2 ) = 0(k 2 ). □ 

Next we recall a result about convergence of step processes towards a diffusion process, see 
Ispany and Pap [29j Corollary 2.2]. This result is used for the proof of convergence (13.71) . 

5.1 Theorem. Let 7 : R + xK4i be a continuous function. Assume that the SDE 

(5.9) dU t = 1 (t,Ut)dW t , teR + , 

has a unique weak solution with initial value Uq = u for all u G R, where (Wt)teR + is a 
standard Wiener process. Let (U t )tm+ be a solution of (15.91) with initial value U = 0. 

For each n G N, let (UJ^)ken be a sequence of random variables adapted to a filtration 
(J r k)k€Z+- Let 

\nt\ 
k=l 

Suppose E ((U£) 2 ) < 00 and E (U% \ ^-1) = f or a ^ n,k e N. Suppose that for each 
T > 0, 



(1) sup 

t€[0,T] 



LntJ 
k=l 



:>>) E E ((t/i") 2 l { |t,;|><,) I 7J_0 A /or oi( 9 > 0, 



fe=l 



P £ 

where — )• denotes convergence in probability. Then U n — > U as n — > oo. 

In fact, this theorem is a corollary of a more general limit theorem, see Ispany and Pap 
Theorem 2.1]. 

Now we recall a version of the continuous mapping theorem. 

For a function / G D(R + ,R) and for a sequence (f n )nen m D(R + ,R), we write 
fn / if (/n)nGN converges to / locally uniformly, i.e., if sup t6[0)T] |/ n (t) - f(t)\ ->■ 
as n — > oo for all T > 0. For measurable mappings $ : D(M + ,M) — > D(R + ,R) and 
$ n : D(R + ,R) -> D(R + ,R), n G N, we will denote by C $i ( $ri ) nsN the set of all functions 

/GC(R + ,R) such that $(/) G C(R+,R) and $„(/„) ^> $(/) whenever f n f with 
/ n G D(R + ,R), n G N. 

For deriving convergence (13. ip from convergence (13. 7p we will need the following version of 
the continuous mapping theorem. 
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5.2 Lemma. Let (Ut)teR+ an d (^4 n )teR + ; fiGN, be stochastic processes with cadlag paths 
such that as n^oo. Let $ : D(R+,R) ->■ D(R+,K) and $ n : D(R+,M) ->■ 

D(R + , R) ; fieN, be measurable mappings such that there exists C C C$,(* n ) neN with C G Poo 
and P(U EC) = l. Then $ n (ZY n ) -A $(W) as n -> c». 

Lemma \5. 21 can be considered as a consequence of Theorem 3.27 in Kallenberg [36J, and we 
note that a proof of this lemma can also be found in Ispany and Pap [291 Lemma 3.1]. 



References 

[1] S. Ahn, L. Gyemin and J. Jongwoo, Analysis of the M/D/l-type queue based on an 
integer- valued autoregressive process. Operations Research Letters 27, 235-241, (2000). 

[2] M. A. Al-Osh and A. A. Alzaid, First order integer-valued autoregressive INAR(l) 
process. Journal of Time Series Analysis 8(3), 261-275, (1987). 

[3] M. A. Al-Osh and A. A. Alzaid, An integer-valued pth-order autoregressive structure 
(LNAR(p)) process. Journal of Applied Probability 27(2), 314-324, (1990). 

[4] P. Belisle, L. Joseph, B. MacGibbon, D. Wolfson and R. du Berger, Change- 
point analysis of neuron spike train data. Biometrics 54, 113-123, (1998). 

[5] P. Billingsley, Convergence of Probability Measures, 2nd ed. Wiley, 1999. 

[6] U. Bockenholt, An INAR(l) negative multinomial regression model for longitudinal 
count data. Psychometrika 64, 53-67, (1999). 

[7] U. Bockenholt, Mixed INAR(l) Poisson regression models: analyzing heterogeneity 
and serial dependencies in longitudinal count data. Journal of Econometrics 89, 317-338, 
(1999). 

[8] K. Brannas and J. Hellstrom, Generalized integer- valued autoregression. Econometric 
Reviews 20, 425-443, (2001). 

[9] K. Brannas and Q. Shahiduzzaman, Integer- valued moving average modelling of the 
number of transactions in stocks. Umea Economic Studies 637, University of Umea, (2004). 

[10] R. A. Brualdi and D. Cvetkovic, A Combinatorial Approach to Matrix Theory and 
its Applications. CRC Press, Boca Raton, FL, 2009. 

[11] A. C. Cameron and P. Trivedi, Regression Analysis of Count Data. Cambridge Uni- 
versity Press, Oxford, 1998. 

[12] M. Cardinal, R. Roy and J. Lambert, On the application of integer- valued time series 
models for the analysis of disease incidence. Statistics in Medicine 18, 2025-2039, (1999). 

[13] N. H. Chan and C. Z. Wei, Limiting distributions of least squares estimates of unstable 
autoregressive processes. The Annals of Statistics 16, 367-401, (1988). 



28 



[14] F. C. DROST, R. V. DEN Akker and B. J. M. WERKER, Local asymptotic normality 
and efficient estimation for INAR(p) models. Journal of Time Series Analysis 29(5), 783- 
801, (2008). 

[15] F. C. Drost, R. v. den Akker and B. J. M. Werker, The asymptotic structure 
of nearly unstable non-negative integer- valued AR(1) models. Bernoulli 15(2), 297-324, 
(2009). 

[16] F. C. Drost, R. v. den Akker and B. J. M. Werker, Efficient estimation of autore- 
gression parameters and innovation distributions for semiparametric integer-valued AR(p) 
models. Journal of the Royal Statistical Society: Series B (Statistical Methodology) 71(2), 
467-485, (2009). 

[17] J. G. Du and Y. Li, The integer valued autoregressive (INAR(p)) model. Journal of Time 
Series Analysis 12(2), 129-142, (1991). 

[18] S. N. Ethier and T. G. Kurtz, Markov Processes. John Wiley & Sons, Inc., New York, 
1986. 

[19] J. Franke and T. Seligmann, Conditional maximum-likelihood estimates for INAR(l) 
processes and their applications to modelling epileptic seizure counts. In: T. Subba Rao 
(Ed.), Developments in time series, pp. 310-330. London: Chapman & Hall, (1993). 

[20] J. Franke and T. Subba Rao, Multivariate first order integer valued autoregressions. 
Technical report. Math. Dep., UMIST, England (1995). 

[21] R. Freeland and B. McCabe, Asymptotic properties of CLS estimators in the Poisson 
AR(1) model. Statistics & Probability Letters 73, 147-153, (2005). 

[22] G. Gauthier and A. Latour, Convergence forte des estimateurs des parametres d'un 
processus GENAR(p). Annales des Sciences Mathematiques du Quebec 18(1), 49-71, 
(1994). 

[23] C. Gourieroux and J. Jasiak, Heterogeneous INAR(l) model with application to car 
insurance. Insurance: Mathematics and Economics 34, 177-192, (2004). 

[24] L. Gyorfi, M. Ispany, G. Pap and K. Varga, Poisson limit of an inhomogeneous 
nearly critical INAR(l) model. Acta Universitatis Szegediensis. Acta Scientiarum Mathe- 
maticarum 73(3-4), 789-815, (2007). 

[25] J. Hellstrom, Unit root testing in integer-valued AR(1) models. Economics Letters 70, 
9-14, (2001). 

[26] R. A. Horn and Ch. R. Johnson, Matrix Analysis. Cambridge University Press, Cam- 
bridge, 1985. 

[27] M. Ispany, G. Pap and M. C. A. van Zuijlen, Asymptotic inference for nearly un- 
stable INAR(l) models. Journal of Applied Probability 40(3), 750-765, (2003). 



29 



[28] M. Ispany, G. Pap and M. C. A. van Zuijlen, Asymptotic behaviour of estimators 
of the parameters of nearly unstable INAR(l) models. Foundations of statistical inference 
(Shoresh, 2000), 195-206, Contributions to Statistics, Physica, Heidelberg, (2003). 

[29] M. Ispany and G. Pap, A note on weak convergence of step processes. Acta Mathematica 
Hungarica 126(4), 381-395, (2010). 

[30] J. Jacod and A. N. Shiryaev, Limit Theorems for Stochastic Processes, 2nd ed. 
Springer-Verlag, Berlin, 2003. 

[31] P. Jeganathan, On the asymptotic behavior of least squares estimators in AR time 
series with roots near the unit circle. Economic Theory 7, 269-306, (1991). 

[32] R. C. Jung, G. Ronning and A. R. Tremayne, Estimation in conditional first order 
autoregression with discrete support. Statistical Papers 46, 195-224, (2005). 

[33] R. C. Jung and A. R. Tremayne, Binomial thinning models for integer time series. 
Statistical Modelling 6, 81-96, (2006). 

[34] M. Kachour, p-order rounded integer- valued autoregressive (RINAR(p)) process. Arxiv: 
0902.1598, (2009). 

[35] M. Kachour and J. F. Yao, First-order rounded integer- valued autoregressive (RI- 
NAR(l)) process. Preprint, University of Rennes 1, (2008). 

[36] O. Kallenberg, Foundations of Modern Probability. Springer, New York, Berlin, Hei- 
delberg, 1997. 

[37] I. Karatzas and S. E. Shreve, Brownian Motion and Stochastic Calculus, 2nd ed. 
Springer, Berlin, 1991. 

[38] A. Latour, The multivariate GINAR(p) process. Advances in Applied Probability 29, 
228-248, (1997). 

[39] A. Latour, Existence and stochastic structure of a non- negative integer- valued autore- 
gressive processes. Journal of Time Series Analysis 19(4), 439-455, (1998). 

[40] S. Ling and W. K. Li, Limiting distributions of maximum likelihood estimators for 
unstable ARMA time series with GARCH errors. The Annals of Statistics 26, 84-125, 
(1998). 

[41] S. Ling and W. K. Li, Asymptotic inference of nonstationary fractional ARIMA models. 
Economic Theory 17, 738-764, (2001). 

[42] I. MacDonald and W. Zucchini, Hidden Markov and Other Models for Discrete-valued 
Time Series. Chapman & Hall, London, 1997. 

[43] T. van der Meer, G. Pap and M. C. A. van Zuijlen, Asymptotic inference for nearly 
unstable AR(p) processes. Economic Theory 15, 184-217, (1999). 



30 



[44] B. P. M. McCabe, G. M. Martin and D. Harris, Optimal probabilistic forecasts for 
counts. Monash University, Working Paper 7/09, 2009. 

[45] E. McKenzie, Some simple models for discrete variate time series. Water Resources 
Bulletin 21, 645-650, (1985). 

[46] E. McKenzie, Discrete variate time series. In: Shanbhag D. N., Rao C. R. (eds) Handbook 
of Statistics. Elsevier Science, 573-606, (2003). 

[47] M. Monteiro, I. Pereira and M. G. Scotto, Optimal alarm systems for count 
processes. Communications in Statistics. Theory and Methods 37, 3054-3076, (2008). 

[48] H. Pavlopoulos and D. Karlis, INAR(l) modeling of overdispersed count series with 
an environmental application. Environmetrics 19, 369-393, (2008). 

[49] P. C. B. Phillips and Z. Xiao, A primer on unit root testing. J. Econom. Surv. 12, 423- 
470, (1998), In: M. McAleer, L. Oxley (Eds.), Practical Issues in Cointegration Analysis, 
Blackwell, Oxford, 1999. 

[50] J. PlCKANDS III and R. Stine, Estimation for an M/G/l queue with incomplete infor- 
mation. Biometrika 84, 295-308, (1997). 

[51] A. M. M. S. Quoreshi, Bivariate time series modelling of financial count data. Com- 
munications in Statistics. Theory and Methods 35, 1343-1358, (2006). 

[52] D. Revuz and M. Yor, Continuous Martingales and Brownian Motion, 3rd ed., corrected 
2nd printing. Springer- Verlag, Berlin, 2001. 

[53] N. Rudholm, Entry and the number of firms in the Swedish pharmaceuticals market. 
Review of Industrial Organization 19, 351-364, (2001). 

[54] M. E. Silva and V. L. Oliveira, Difference equations for the higher-order moments 
and cumulants of the INAR(l) model. Journal of Time Series Analysis 25(3), 317-333, 
(2004). 

[55] M. E. Silva and V. L. Oliveira, Difference equations for the higher-order moments and 
cumulants of the INAR(p) model. Journal of Time Series Analysis 26(1), 17-36, (2005). 

[56] I. Silva and M. E. Silva, Asymptotic distribution of the Yule- Walker estimator for 
INAR(p) processes. Statistics & Probability Letters 76(15), 1655-1663, (2006). 

[57] F. Steutel and K. van Harn, Discrete analogues of self-decomposability and stability. 
The Annals of Probability 7, 893-899, (1979). 

[58] F. W. Steutel and K. van Harn, Infinite Divisibility of Probability Distributions on 
the Real Line. Dekker, New York, 2004. 

[59] P. Thyregod, J. Carstensen, H. Madsen and K. Arnbjerg-Nielsen, Integer val- 
ued autoregressive models for tipping bucket rainfall measurements. Environmetrics 10, 
395-411, (1999). 



31 



[60] C. H. WeiB, Thinning operations for modelling time series of counts — a survey. Advances 
in Statistical Analysis 92(3), 319-341, (2008). 

[61] J. Zhou and I. V. Basawa, Least-squared estimation for bifurcation autoregressive pro- 
cesses. Statistics & Probability Letters 74, 77-88, (2005). 



32 



